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5-1 ' Abstract 

, The Rastall's theory is a modification of the General Relativity theory leading to a different expression 

\^ . for the conservation law in the matter sector compared with the usual one. It has been argued recently that 
such a theory may have applications to the dark energy problem, since a pressureless fluid may lead to an 

■ accelerated universe, hi the present work we confront the Rastall's theory with the power spectrum data. 

\ The results indicate a configuration that essentially reduces the Rastall's theory to General Relativity, unless 

. the non-usual conservation law refers to a scalar field, situation where other configurations are eventually 

I ' possible. 
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I. INTRODUCTION 



A large number of cosmological observational data requires the existence of two exotic compo- 
nents in the matter content of the universe, dark matter and dark energy l|. Both constitute the 
so-called dark sector of the cosmic budget. Dark matter is necessary, for example, to explain the 
formation of structures in the expanding universe. Dark energy is necessary, on the other hand 
to explain the present stage of the accelerated expansion of the universe. Dark energy must have 
a negative pressure in order to induce the acceleration of the universe. This is a quite strange 
property, and one of the most important question today in theoretical physics concerns the nature 
of such an exotic fluid. All candidates to describe the dark energy component face serious problems 
and drawback. 

Another possibility to explain the observational data is to consider that General Relativity is not 
the true gravitational theory. Some modifications of General Relativity may lead to an accelerated 
universe even if only usual types of matter are taked into account. A very fashion proposal is the 
f{R) theories, a non-linear generalization of the Einstein-Hilbert action . For a review of the f{R) 
theories and their present status, see reference Q]. In fact, the unusual properties of dark energy 
motivate the search of many other alternatives to explain the observational data. 

Recently, it has been investigated if the Rastall's theory may be a viable alternative to the intro- 
duction of dark energy in the General Relativity context 3]- The Rastall's theory is a modification 
of the General Relativity, proposed in 1972 jj]. It implies a change of the Einstein's equation that 
mounts out to a modification of the conservation law for the energy-momentum tensor. One of 
the motivations to this proposal is the fact that the usual conservation law is only firmly tested in 
special relativity. The theory contains a free parameter A such that A = 1 implies General Relativ- 
ity. In this way, it can be considered as a deformation of General Relativity. The modification in 
the conservation law can lead to new effects compared with General Relativity. Depending on the 
value of the parameter A, a pressureless matter, for example, can induce an accelerated expansion. 

An important drawback of Rastall's theory is the absence of the Langrangian formulation. But, 



it is possible that it can obtained from an action principle using Weyl's geometry 
argued in reference 



It has been 



that the Rastall's theory may appear even in the context of Riemannian 
geometries by a redefinition of the energy-momentum tensor. In any case, the Rastall's theory 
implies that a given fiuid, with a specific equation of state, may have a very different effective 
equation of state when interpreted in the context of General Relativity. This properties has called 
the attention to this proposal as an alternative to dark energy. 
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In reference [3|, the Rastall's theory has been confronted agains the supernova type la data. 
The result indicates that it is competitive with respect to the ACDM model, but the mass density 
is very high, around ~ 0.5. This suggests that the theory may face troubles when tested against 
the mass agglomeration phenomena. 

In the present work, we intend to test the Rastall's theory using the power spectrum observa- 
tional data. This implies a perturbative study. A perturbative study of the Rastall's theory has 
already been performed in reference [Ti]. The main conclusion was that the theory behaves also 
at perturbative level like General Relativity but with an effective equation of state. In particular, 
if the fluid has an effective equation of state necessary to induce the accelerated expansion of the 
universe, it will present a negative squared effective sound speed. Hence, the resulting scenario is 
plagued with instabilities. Moreover, in order to compare the theoretical predictions of the Rastall's 
theory with the power spectrum data, it seems unavoidable to introduce a two fluid model. This 
represents a problem in the original Rastall's theory. A possible way out is to consider a fluid which 
obeys the conservation law of Rastall's theory, with another one that conserves separately in the 
usual way. This may be justified if the former fluid is in fact a field, like a scalar field, that obeys 
a modified Klein-Gordon equation. In fact, if both fiuids have a hydro dynamical representation, 
the power spectrum predicts A = 1 as it will be described later in this work - that is, the theory 
reduces to General Relativity. If the Rastall's fiuid is a scalar field, the situation is more complex, 
and even if A = 1 remains favored, other possibility appears. 

This paper is organized as follows. In next section, we present the Rastall's theory, deriving 
some cosmological relations. In section III, the matter power spectrum is determined in the hy- 
drodynamical representation. In section IV, the matter power spectrum is determined for the case 
where one of the fiuids is representing by a self-interacting scalar field. In section V we present our 
conclusions. 



II. THE FIELD EQUATIONS 

Originally, the fundamental equations of the Rastall's theory were written as 

R/iu ~ '^dfJ.v^ — I^Tf^u, (1) 

j.f,u ^ }—Ar;'^, (2) 
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The General Relativity theory, with the usual conservation of the energy-momentum tensor, is 
re-obtained when A = 1. These equations may be re- written as 

Rfiu - -j^g^iuR = ^^^^Jiu — '^^—^ — ^Ati/^j, (3) 
T'^^/. = (4) 

where 

3A-2 

Again, 7 = 1 (corresponding to A = 1) implies General Relativity. Equivalently, one can recast 
these equations under the following form: 

RtJiv = f^^tJiu -'^^-Y-g^lvT^, (6) 

T'^^/. = (7) 

In analysing the perturbed Rastall model, we must consider a two fluid model. The reason 
is that baryons clearly exists, and there are good observational evidences that baryons can be 
modelled by a pressureless matter with an approximately zero sound velocity. This indicates that 
baryons can not be introduced directly in the framework of equations (|l|2p by simply adding a 
new energy-momentum tensor in the rhs of ([1]): if the original framework of the Rastall's theory 
is preserved, even if we set a fluid with zero pressure the resulting sound velocity is not zero. 
Alternatively, a conserved baryonic energy-momentum tensor could be added to the rhs of ([1]), 
but this would imply a presence of an inhomogeneous term for the exotic fluid obeying ([2|). This 
inhomogeneous term leads to a negative energy component of the exotic fluid, which dominates 
either in the past or in the future. For this reason, we will consider equations (j3|4p the fundamental 
framework of the two fluid model including baryons; it implies the addition of a second energy- 
momentum tensor which obeys the usual conservation law. This additional energy-momentum 
tensor will represent the baryonic component. 

Hence, we will consider a cosmological model with two fluids, one obeying the Rastall's frame- 
work, with no usual conservation of the corresponding energy-momentum tensor, and the other 
obeying the traditional conservation law of general relativity. Under these conditions, the field 
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equations read, 

R,u = ^{t^. - '^-=^9,uT^Y^TjZ, (8) 

TH",^ = ^Tf, (9) 

Tr^:,, = 0. (10) 

The subscripts (superscripts) x and m designate the exotic and matter (baryonic) components, 
respectively. 

The universe is homogenous and isotropic at least at scales greater than 100 Mpc. Hence, 
at sufficiently large scales, it can be represented by the Friedmann-Lemaitre-Robertson- Walker 
(FLRW) metric, 

ds' = - a(t)2|-^!^ + r^{de^ + sin2 0d,/.2)|, (11) 

where a{t) is the scale factor and k is the curvature of the spatial section. Introducing this metric 
in the field equations (j8l9|10p . and specializing for the flat case {k = 0), we obtain the following 
equations of motion: 

-j = g S 3 - 7 - 3(1 - 7)tJx ^Pi^ + -Pm, (12) 

^ = -^|7 + 3(2 -7)wx|px - ^/Jm, (13) 
(X 'y — 1 

Px + 3-(l + UJx)Px = ~^(1 ~ 'i^x)Px-, (14) 

d 

Pm + 3-Pm = 0. (15) 

a 

The two last equation can be integrated leading to, 

Px = Pxoa^-'+^(-''^^-- , (16) 

Pm = PmOa~^. (17) 

Remark that all this formulation is equivalent to a traditional General Relativity, where the perfect 

fluid would have an effective equation of state given by, 

_ 7 - 1 + (5 - 87)0;^ 
"^^^ - 3-7 + 3(7-1)0;.- ^''^ 

It is interesting that = — 1 implies uJeff = the cosmological constant case is a fix point 

in this stuff. The reason for this fixed point is easily understood inspecting ([2]): it corresponds 

to a de Sitter (or anti-de Sitter) space-time, whose curvature is constant; consequently, the usual 

conservation law is recovered. 



5 



Using the modified the flat Friedmann's equation, and defining, as usual. 



the following equation must be obeyed today: 



^xo - -^jp , ^^0 - -^jp , (19) 



l = ^|^|3_^_3(l_^)a;,|+[7o. (20) 

Remark that to define the density parameters the newtonian cosmological constant is employed. 
It could be used an effective cosmological constant, leading to a numerical different value for the 
density parameter, but without changing the general framework Q]. 

III. PERTURBATIVE ANALYSIS: THE FLUID DESCRIPTION 

Let us consider in Rastall's theory the behaviour of a given fluid characterized, for example, by 
an equation of state p = cop, with oj constant. The predictions of the Rastall's theory, in one fluid 
description, is equivalent for the background point of view, to the predictions of General Relativity 
when a fluid with an equation of state of the type p = tJeffP, ^ and Wg// being connected by the 
relation (llSp . In this sense, with this identification, one theory can be mapped into the other. 

At perturbative level, however, this equivalence is not evident. In reference (Ql)) a perturbative 
study has been carried out, considering just one fluid, obeying the framework of the Rastall's 
theory. The fluid description was kept all along the calculation. The final equations reveal the 
same results of General Relativity, provided that the identification (jlSp is made. The equivalence 
remains at perturbative level as far as the fiuid description is used. 

One of the most important problems today in cosmology is the description of dark energy. 
Dark energy requires negative pressure and, at perturbative level, a fiuid with negative pressure is 
unstable at small scales Q]. The Rastall's theory opens new possibility, but as far as we remain 
at the level of a fiuid description, and with a one-fiuid model, it seems that the same problems 
remain: a negative effective pressure is required at background level, and this leads to instabilities 
at perturbative level using the results of reference [71]. 

But, as already briefiy discussed in the Introduction and in the previous section, we can go one 
step further by considering the two fiuid model, one of them representing the baryons. This allows 
to consider the observational data to restrict the model. In fact, the power spectrum concerns the 
baryonic component. In other words, it concerns a fluid with zero effective pressure, what assures 
the gravitational collapse, leading to the formation of local structures. The other fluid will follow 
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the Rastall's framework. Its considered as the dark component of the universe. To be specific, this 
exotic fluid wiU be taken having zero pressure, leading to an effective equation of state parameter 

^eff = (21) 

In this situation, acceleration of the universe can be achieved if 7 < or 7 > 3. Moreover, relation 
(pO|) reads now, 

1 = ^^x-o^ + ^0. (22) 

Let us now consider the perturbation of the two fluid models. We will work in the synchronous 
gauge. Since the modes of interest for the matter power spectrum are well inside the horizon, the 
results are not sensitive to the use of one gauge or another, or even a gauge invariant formalism. 

In the synchronous gauge condition, we introduce the perturbations in the metric and matter 
functions, 

Px = PxO + SPx , Pm = PmO + ^Pm , Um = UmO + , (23) 

Px = PxO + Spx , Pm = PmO + ^Pm , Ux = UxO + , (24) 

9fii^ = gfiu + hf,u , h^o = 0. (25) 

In expressions (|23ti25p . the sub(super)scripts "0" indicate the background functions, and 6pm, ^Px, 
5pm, Spx, 6um, Sux, indicate the perturbed quantities in density, pressure, four- velocity and 
metric, while /i^o = defines the coordinate condition. Since both fluids have zero pressure, we 
fix Px = Pm = 0. Moreover, since we will not taken into account entropic perturbations (which 
nevertheless may lead to new interesting effects in the present framework) we can fix 5px = 5pm = 0. 

Introducing these perturbations and remaining at the linear level, we obtain the following 
coupled equations for the perturbed quantities: 

Sm + 2-6m = ^jil.xoa~^Sx + ^Qoa~^5m, (26) 
2 

^x = (Ox-Sm), (27) 

3-7 



9-57d _ 1 - 7 ,,2 

"x — 7: 1^ —n: 

'J a 2, a'^ 



+ ^^-Jx = if ' (28) 

(29) 



h 
2" 



In these expressions, we have the following definitions: 



^ , 5m=^-^ , e = d,6u, , h = ^i^ . (30) 



In order to integrate numerically perturbed equations (j26ll28p it is more convenient to use the 
scale factor a and as the dynamical variable, since it is directly connected with the dimensionless 
redshift quantity through z = — 1 + ^ (fixing today ao = 1). In terms of this new variable, we 
obtain the following equations: 

A" f'(a)\ , 3 _ 6 3 no _3 

= -3^(:^-^0' ^^^^ 



9-576*^ _ 1-7A;2 
3-7 "a ~ 2~fcf a2/(a)' 



''^ + ^^ = -^T27jfrr., (33) 



(34) 



where ko is the wavenumber associated to the Bubble's radius, and it has been defined the function, 



fia) 



3-7 ^ , f^o 



1/2 



2 a 

To obtain a prediction, we compare the matter power spectrum, defined by. 



(35) 



V = 5i, (36) 
n 

with the observational data of the 2dFGRS observational program 19]. To fix the initial condition 

nn 

we use the BBKS transfer function [10|, [11|], and apply the prescription described in the reference 
[l^ . We use the statistical parameter defined by 

X =2^ -2 > (37) 

where "P*^ is the ith observational data, af being the error bar, is the corresponding theoretical 
prediction. The probability distribution function (PDF) is given by 

P(7,17o) = ^e-^'/2^ (38) 

where ^ is a normalization constant. As indicated, the PDF depends on two free parameters, the 
matter density f^o and the 7 parameter which characterizes the deviation from the Einstein theory. 
Marginalizing (integrating) in one of the parameters, we obtain the one-dimensional PDF for the 
remaining free parameter. The results are shown in figure 1. The probability is highly concentrated 
around 7 = 1, which corresponds to the Einstein theory - it admits a slight deviation if the matter 
parameter is high, almost without the exotic fluid. This result can be easily understand: as far as 
the parameter 7 differs from the General Relativity value 7 = 1, oscillations or even exponential 
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FIG. 1: On the left, the two-dimensional probability distribution function (PDF) for the matter density flo and the 
7 parameter. In the center the one-dimensional PDF for the 7 parameter, and on the right the one-dimensional PDF 
for the matter density parameter r^o. 

behaviour (depending if 7 is lesser or greater than 1, respectively) is induced in the exotic fluid, 
and this is highly transfered to the matter fluid. Such behaviour is not observed in the matter 
spectrum. Remark that if ^IxO ~ 0, implying ~ 1, larger deviations from 7 = 1 are possible. 

Hence, we can conclude that the value 7 = 1 is the prediction obtained from the matter power 
spectrum, with a precision of the order of 10~^ (see the graphics). In fact, the best fitting is 
achieved by 7 = 1 and ~ 0.79, = 0.38 per degree of freedom. The best fit model ACDM 
model has the same P^r degree of freedom. 



IV. PERTURBATIVE ANALYSIS: THE SCALAR FIELD DESCRIPTION 

The result described above is restricted to a fluid formulation of both components, that repre- 
senting matter and that representing the exotic fluid. This exotic fluid must, more precisely, be 
represented by a field than a fluid - it is very difficult that a usual fluid could present the exotic 
behaviours connected with dark matter and dark energy, or even follow the new conservation law 
dictated by the Rastall's theory. In this sense, let us consider the most simple field description in 
cosmology, that of a self-interacting scalar field. The energy momentum tensor of this field is given 
by, 

Tf^u = 4';fi4';u - ]j^g^lv4>'''(t>■,p + giivV{(l)). (39) 



9 



If we impose the usual conservation law for the energy-momentum tensor T^'^-^ = 0, we obtain the 
usual Klein- Gordon equation: 

VpV^ = -V^{^). (40) 

However, if now the conservation law must read as in equation ([7]), the scalar field must now obey 
the following equation: 

V,VPct> + (3 - 27)F^ = (1 - 7)^^^. (41) 
In the two fluid model, the corresponding "Einstein's" equation takes the form, 

From the scalar field description it comes out some similarities of the Rastall's theory with the K- 
Essency models [l4]. K-Essence models may also be plagued with negative sound velocity problems 
what strength the mentioned similarities [isf. 

One clear advantage in using the scalar field representation is that it usually avoids the problems 
of the fluid representation of components with negative pressure. However, while this is clear in 
the usual case, with the Einstein's equation coupled to the ordinary Klein-Gordon's equation, this 
is less clear in the framework of the Rastall's theory. In order to investigate the constraints from 
mass power spectrum for this field formulation of the Rastall's theory, we take the coupled system 
matter -|-scalar field-t-Rastall's gravity, written in a convenient way for the perturbative analysis in 
the synchronous coordinate condition: 

Rp.v = ^■,^<t>;u + ^-^g^u4>-A''' - (3 - '2l)9^uV + 87tg(^T^, - Igt^uT^ , (43) 

V,V^<^ = -(3 - 2j)V^ + (1 - i f'fJX'^ (44) 
T^".^ = 0. (45) 



Now, let us suppose that the scalar field has zero pressure. Since the energy and the pressure 
of the scalar field is given by 

= y + ^('/'), (46) 
P<A = f-n0), (47) 

the condition of zero pressure implies (j)'^/2 = V, leading to = <p'^. Prom the point of view of the 
background, there is no difference between the fluid and scalar field approach. But at perturbative 
level the difference is significative. 
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In the synchronous coordinate condition, the perturbed equations corresponding to the system 
described by (|43II45|) . are the following: 

6 + 2-6-^^5 = {3- 7)<A^ - (3 - 27)^^^, (48) 



^^ + 3_^+|_ + (3_27)y^^j^' = </,5, (49) 

where = 5(j) and 5 is the density contrast of the matter component. Using now the scale factor 
as the variable, the above system of equations take the following form: 

2 , /'(a)l i 3 !!„ , ,„ ,„ Vt 



where 



-3 



^ = y^oj^, (52) 
V = ^n,,oa^, (53) 

= f > (54) 
V,, = '-{V,. (55) 

The function /(a) is defined as before. 

We perform the same statistical analysis as before, again imposing the initial conditions using 
the BBKS transfer function. The results are shown in figure 2. The main difference is that there 
is now two relevant regions in the space parameter: one around 7 = 1, with a low density, and 
the other near 7 = 0, but positive, extending from the low to high density. The region around 
7 = 1 has high probabilities, but it is smaller; the region near 7 = has lower probabilities, but 
extend to a large region. The consequence is that, after marginalization, there is two peaks in the 
one-dimensional probability for 7: one near 7 = 1 and another near 7 = 0. The second peak is 
higher. We think that this is an effect of the larger probability region around 7 = 0, which seems to 
compensate the higher probabilities around 7 = 1. There are also two peaks in the one-dimensional 
PDF for JIq) one near JIq = and another near = 1- The best fitting is achieved for 7 = 1.02 
and r^o = 0.72 with a ^ per degree of freedom equal 0.30 better than the ACDM model. The 
PDF goes quickly to zero for 7 < 0. Still concerning the peak near 7 = 0, the parameters for 7 
and r^o in this region implies a around 0.33, high compared with the for 7 ~ 1, but smaller 
than the ACDM best fitting. We may ask about the statistical relevance of this second peak. 
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FIG. 2; On the left, the two-dimensional probability distribution function (PDF) for the matter density CIq and the 
7 parameter. In the center the one-dimensional PDF for the 7 parameter, and on the right the one-dimensional PDF 
for the matter density parameter Qq. 



V. CONCLUSIONS 

In this work we have investigate the Rastall's proposal of modification of General Relativity 
with respect to the problem of structure formation. This proposal is equivalent to change the 
usual conservation law of General Relativity. In some sense, it can be seen as a modification of 
the equation of state of a given fluid from the dynamical point of view. In this sense, the Rastall's 
theory can be interesting in order to obtain an accelerated expansion of the universe without 
introducing exotic fluids. For example fluids with positive or null pressure may induce a dynamics 
typical of fluids of negative pressure. This behaviour has already been remarked, for example, in 
the Brans-Dicke theory 15|. 

Our results indicate that the Rastall's theory faces many problems at perturbative level. Con- 
sidering a homogenous and isotropic universe for example, the effective equation of state of the 
background is entirely reproduced at perturbative level, leading to high instabilities when this ef- 
fective equation of state implies negative pressure. One way out is to consider a two- fluid model, 
one of the fluids obeying the usual conservation law and the other one following the Rastall's 
prescription. The results point out to a configuration which still reduces the Rastall's theory to 
General Relativity. 

When we keep the framework of two components which obeys different conservation laws, but 
with a self- interacting scalar field playing the role of the "Rastall's fluid" , more interesting features 
appear. In particular, the Rastall's proposal can be seen as a modification of the Klein-Gordon 
equation, similarly to what happens in the K-Essence theories [l^. Again the configuration cor- 
responding to General Relativity is favored, but other configurations with 7 ~ are also possible 
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even if some statistical subtleties appear. All these considerations seem to indicate to the specificity 
of General Relativity and rule out the Rastall's theory. 

A possible way out to save the Rastall's proposal is to consider the modification of the usual 
conservation law as a manifestation of quantum effects (like particle production) in the spirit of 
reference The particle production in an expanding universe may lead to new terms in the usual 
conservation law. Such possibility remains to explore. In any case, it seems clear that structure 
formation asks for a fluid we behaves in the background and at perturbative level with zero effective 
pressure (otherwise no mass agglomeration can effectively occur) , and this poses a serious problem 
in the original framework of the Rastall's theory. In some sense, this has already been remarked 
in reference Q], forcing the authors of that work to use just one fluid to fit the supernova type la 
data. But, such procedure seems to be impossible concerning the structure formation problem. 

Other observational methods, like CMB, may be used to constrain better the Rastall's theory. 
But, the results obtained in the present work seem to us to be strong enough to point out the 
difficulties that the Rastall's theory face when confronted with observational data. 
Acknowledgements: We thank Oliver Piattela, Winfried Zimdahl and Hermano Velten for their 
remarks and criticisms on the text. We thank also CNPq (Brazil) and FAPES (Brazil) for partial 
financial support. 



[1] R.R. Caldwell and M. Kamionkowski, Ann. Rev. Nucl. Part. Sci. 59, 397(2009). 

[2] A. de Felice and S. Tsujikawa, f(R) theories, arXiv:1002.492l 

[3] M. Capone, V.F. Cardone and M.L. Ruggiero, Accelerating cosmology in Rastall's theory, 

larXiv:0906.413 9( 

[4] P. Rastall, Phys. Rev. D6, 3357 (1972). 

[5] L.L. Samley, Class. Quant. Grav. 10, 1179 (1993). 

[6] L. Lindblom and W.A. Hiscock, J. Phys. A, 1827(1982). 

[7] J.C. Fabris, R. Kerner and J. Tossa, Int. J. Mod. Phys. D9, 111(2000). 

[8] J.C. Fabris and J. Martin, Phys. Rev.D55, 5205(1997). 

[9] S. Cole et al, Men. Not. Roy. Astron. Soc. 362, 505(2005). 

[10] N. Sugiyama, Astrophys. J.Suppl. 100, 281(1995). 

[11] J.M. Bardeen, J.R. Bond, N. Kaiser and A.S. Szalay, Astrophys. J. 304, 15(1986). 

[12] J.C. Fabris, I. Shapiro and J. Sola, JCAP 0702, 016(2007) 

[13] R. de Putter and E.V. Linder, Astropart. Phys. 28, 263(2007) 

[14] C. Armendariz-Picon, V.F. Mukhanov and P.J. Steinhardt, Phys. Rev. Lett. 85, 4438(2000). 



13 



[15] A.B. Batista, J.C. Fabris and R. de S Ribeiro, Gen. Rel. Grav. 33, 1237(2001). 
[16] S.H. Pereira, C.H.G. Bessa and J.A.S. Lima, Quantized fields and gravitational particle creation 
expanding universes, larXiv:091 1 .06221 



14 



